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Abstract 

A simple exact-exchange density-functional method for a quasi-two- 
dimensional electron gas with variable density is presented. An analytical 
expression for the exact-exchange potential with only one occupied subband 
is provided, without approximations. When more subbands are occupied the 
exact-exchange potential is obtained numerically. The theory shows that, 
in contradiction with LDA, the exact-exchange potential exhibits disconti- 
nuities and the system suffers a zero-temperature first-order transition each 
time a subband is occupied. Results suggesting that the translational sym- 
metry might be spontaneously broken at zero temperature are presented. An 
extension of the theory to finite temperatures allows to describe a drop in 
the intersubband spacing in good quantitative agreement with recent experi- 
ments. 



1 



Density-functional theory (DFT) [T[| has been one of the most successful approaches to 
the problem of interacting particles. The use of the local-density-approximation (LDA) 
of DFT and related approximations has been widely applied to the study of atoms, 
molecules, clusters, nuclei, and condensed matter systems. |||]] At the same time, steady 
improvements on the experimental techniques have made possible the production of a great 
number of semiconductor artificial structures that confine the 3D electron gas in one or 
more dimensions. In particular GaAs/ALjGai-^As heterostructures have attracted intense 
attention, because they are particularly well grown, their conduction band is simple to 
describe, and they allow the continuous change of the electron density. Because typical 3D 
electron densities in the quantum well (QW) are very small, the conduction band dispersion 
can be approximated, without significant error, by a second order expansion around the 
minimum (effective- mass approximation). || As a consequence, the effective electron gas 
formed in the nanostructure can be treated in the frame of Kohn-Sham DFT. All the 
theoretical techniques developed around DFT and its approximations, such as the LDA, time 
dependent LDA (TDLDA), spin-density LDA (LSDA) can be tested in these structures which 
are much easier to modify and control than natural materials ||. Therefore, semiconductor 
nanostructures offer the possibility to test over a broad range of conditions the validity of 
the theoretical tools that are used to describe matter in general. 

In this paper we present an approach that goes beyond LDA: an exact-exchange ||[7] 
density-functional theory for the case of a quasi- two-dimensional electron gas (2DEG). Be- 
cause we study a metallic system with several subbands occupied, our theory is different 
than previous exact-exchange (XX) procedures aimed to the study of atoms or bulk semi- 
conductors J7|. Moreover, the number of particles is not fixed but it is allowed to fluctuate 
||. When a single subband is occupied we provide an analytical expression for the XX 
potential without approximations. When more subbands are occupied the XX potential is 
obtained up to numerical precision. The implementation of the method is simple enough to 
replace the state of the art method for 2DEG which is currently based in the LDA. || We 
prove that the XX potential is, in general, a discontinuous function of the density which can 
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not be accounted by LDA. Therefore, the XX approach yields qualitatively new physics: at 
zero temperature a first-order transition occurs every time a new subband is occupied. In 
the transition region, the theory suggest a new phase in which the translational symmetry 
of the system might be broken. An extension of the theory to finite temperatures allows us 
to describe a drop in the intersubband spacing when the first-subband is occupied in good 
quantitative agreement with recent experiments ||. 

This theory can be applied to any system with translational symmetry in a plane. For 
instance, a QW grown epitaxially as shown schematically in Fig. |l]a. In these systems, it is 
possible to confine an electron gas changing the semiconductor in the growth direction z. If 
the larger gap semiconductor is doped with donors, it provides electrons to the trap formed 
by the smaller gap semiconductor, until the charge-transfer field equilibrates the chemical 
potentials in the QW and the doped region (see Fig. |I|b ). The charge transferred from 
the doped reservoir to the QW can be tuned by an external electric field 0. Assuming 
translational symmetry of the 2DEG in the x — y plane (area A) , and proposing accordingly 
a solution of the type (j)^ (r) = exp(z k ■ p)£? (z) / \fA the ground-state-electron density can 
be obtained by solving a set of effective one-dimensional Kohn-Sham equations of the 
form 

a (*) = efff (z) , (1) 

where h 2 /(m*a,Q 2 ) = 2 Ry* and are the units of energy and length, being m* the electron 
effective mass and Oq the effective Bohr radius. £f (z) is the wave function corresponding to 
an electron with a spin projection a (| or j), subband index i, and eigenvalue e° . In Eq.(P 
Vks(z, a) = V ex (z) + Vh(z) + Vxx(z , a)+V c (z, a) is the spin-dependent Kohn-Sham potential 
given as sum of the external, Hartree, local XX, and correlation potentials, respectively. 
The external potential is given by the sum of the epitaxial potential plus any external 
electric field. The zero-temperature 3D electron density n(z) = Y,cr n (. z i°') — J2 £ rT <^f 1 ~ 
£ i) \£i ( z )\ 2 l^i n ( z i a ) being the fraction of a polarized electrons, /i is the chemical potential 
(or Fermi level) of the system, which is determined by the electrostatic and thermodynamic 
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equilibrium with a reservoir. Vh{z) is the solution of the Poisson equation. We approximate 
the correlation potential V c (z, a) as a function of the local spin density ]TU|. It remains to 
define the XX potential Vxx- 

The exchange energy Ex [n(z, a)] of a 2DEG can be obtained through a Slater determi- 
nant constructed with the (occupied) self-consistent solutions of Eq. (1), as follows 

£x = -E dz dz' tf(z, z')^(z', z)F£(z, z'\ (2) 

where tf(z,J) = F$»{z, z') = £ S ^-/^/}^ a p) , # = - 

ef) (/i — ^f) 1 ^ 2 , is the step function, and J n (x) stands for the cylinder Bessel function 
of order n. The XX potential Vxx (z, a) can be obtained from 

5E x a f j ,Wks(z',o-) 



V X x(z,a) = X = I dz 



,. ., Sn (z,a) 



x 



dz 



6E X 5 ^\ z 



+ c.c. 



■ SE X 6k$ l 
+ 5kf5V KS (z',a) < 



5tf(z")5V KS (z>,a) 

The first term in the r.h.s. of Eq.(|3p comes from functional derivatives with respect to 
the "shape" of the wave function j7j, while the second term is a result of changes in the 
occupation of the subbands and has not been considered in previous XX treatments for fixed 
particle-number systems. SVks (z 1 ,a) /5n(z,a) = x~ x {z,z'} is the inverse of the operator 
Xa(z, z') = Sn (z,a) jSV KS {z') given by 

? 

Alt A ir % K~,~j-»y.~i~, . —.j ^ 



Xa(z, Z>) = £ { [tf (*, z')G ia (z, Z') + C.C.] - ^ \^( Z , z')f \ , (4) 



being G ia (z, z') = Y.^i)^{z\z)/{e1 - e°), and /(ef) = 1/ [1 + exp((ef - »)/k B T))] the 
Fermi occupation factor [111]. The first term in Eq.(EI) comes from first-order perturbation 



theory ||12|| , whereas the second term results from first-order perturbation theory and the 
thermodynamic equilibrium between the 2DEG and the reservoir that fixes a common chem- 
ical potential /i allowing the change of the number of particles. Indeed, without this second 
term, the operator x<r{z> z> ) in general cannot be inverted because it is singular Eq.@ 



could be brought to an alternative expression, which allows the discussion of our approach 
in the context of the Optimized Potential Method (OPM). |jTB[ For this, we multiply Eq.(|3]) 
by Xcr{ z i z ')i an d integrate over all z. Proceeding this way, it could be rewritten as 

J» CT ) 2 J dz'{[V xx (z', a) - v X4 (z', a)] G ia {z', z)tf{zf, z) (5) 

2tt 8 E x ' 



+cc}-2j:ierwr 



VxxAa) 



Akf dkf 



0. 



Here, vx,i(z,a) = [4tt / 'A(kp) 2 ^(z)*]SE x / 'S^(z), and V XX) i{a) * s the diagonal matrix ele- 
ment of Vxx(z,cr) with £f(z). This integral equation for Vxx(z,a) is the generalization of 
the OPM to our open configuration, where the system is free to exchange particles with the 
surroundings; for a closed system (fixed number of particles), the last term on the r.h.s. of 
Eq.(||) is zero, and the integral equation reduces to the one of the standard OPM for atoms 



and molecules. [13] Some consequences of Eq.@ are worth of address: a) the solution for 
Vxx( z , &) is univocally determined, including the (possible) presence of a constant shift C; b) 
its solution in the one-subband case i = is immediate: Replacing V xx (z, a) = vx,o(z, cr)+C 
in Eq.(^) and solving for C, we found 



F^(z u z 2 ) F^( Zl ,z 2 ) 



V xx (z, °) = —^y dz^Uzi, zi)F£(zi, z) + k$ J d Zl dz 2 \<pl(z x , z 2 ) 

(6) 

Using Eq.(||) it can be shown that JV xx (z, a) (z)] 2 = dE x /dN a = /i x , with E x given 
by Eq.(|), N a = J dr n(z, a) the total number of electrons with spin a, and [f x the exchange 
contribution to the total chemical potential \i. Accordingly, the Janak theorem |L4j] is explic- 
itly satisfied in our XX-DFT formalism. Besides, when Eq. @ is particularized to the strict 
2D limit by imposing the condition [£q (z)} 2 — > 5{z) [8{z) being the Dirac delta function], we 
obtain V xx (0, a) = — 2k Q p /it, which is exactly the well known value of the exchange contri- 
bution to the chemical potential for an homogeneous 2D electron gas. In other words, Eq. (0) 
contains the exact 2D limit of the exchange potential. JT3[ Finally c) ; for two (or more) occu- 
pied subbands, the analytical exact solution of our OPM integral equation is unknown. Some 
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progress can be achieved following the Sharp-Horton (SH) or KLI approximation [13| of the 



orbital Green function Gi a (z,z'). It amounts to replace the denominators (ef — ej) by an 
orbital independent average Ae a ; substituting G ir7 (z, z') ~ 5(z — z') — Yj% Wl{ z i z ')\ 2 /Ae a 
in Eq.(^]) one obtains an explicit solution for Vxx(z,<r), which will be given elsewhere. |I6(| 
For the one-subband case, the KLI approximation leads to the exact result, given by Eq.([]). 
But for the multi-subband case {i ^ 1), and because of the second term in Eq.(|), the 
explicit solution for Vxx{z,a) becomes Ae a dependent, in contrast with the situation for 



atoms and molecules. [13 



As proceeding along this line of work would had forced us to introduce a new and 
unknown scale of energy (Ae^), we studied the multi-subband case by directly solving 
Eq.(4). The fundamental ingredient for this direct approach is an analytical limit for 
Xo- 1 {z, z'). For T = and in the case where only the ground state subband is occu- 
pied Xo*( z , z ') = Ei dfipf (z) {z') ,where <% = -2 , df >0 = -(ef - e^)/(ji - ej), and 

(z) = £f (z) /^Q (z) . Replacing this and Eq.(0) into (0) we get again, by a different 
method, Eq.(||). For many occupied subbands, one can obtain Vxx(z,a) using Eq.@ be- 
cause x& l { z i z ') can be obtained recursively in terms of Xoa ( z ? z ') f° r an y number of bands 
and temperature, by using the Sherman-Morrison method. [|T7|] Thus for more than one sub- 
band we have evaluated Vxx{z, o~) up to numerical precision. In order to test the numerical 
method we verified that, in double quantum well systems with two occupied subbands, the 
single subband analytical limit for the XX potential in one well (Eq.(§)), is nicely reproduced 
numerically as we increase the barrier between the wells. 

Fig.^b gives the full self-consistent potential Vks(z) and the squared wave functions 
[£? ( z )] 2 corresponding to the first four states. The results were obtained for a QW of GaAs 
in Al x Gai_ x As with a band offset of 220 meV and a width of 245 A. The structure is doped 
on one side at 135 A from the left QW edge with a 3D dopant density of 2 x 10 18 /cm~ 3 . 
Without external electric field and at zero temperature the ground state is paramagnetic, 
the intersubband spacings are E m = — Eq = 26.78 meV, E 02 = 60.58 meV, E 03 = 108.36 
meV, while \i — = 24.14 meV is just below — e$ and only one subband is occupied. 

In Fig.|| we compare the potential V X x( z i a ) corresponding to Fig.|l]b [obtained using 
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Eq.(^)] with the exchange potential in LDA V^ DA (z, a) = — [6n(z,a)/ir] for the same 
XX density. Note that although V^ DA and V xx have a similar amplitude in the QW 
region, their respective asymptotic behaviors are completely different: while V^ DA goes 
exponentially to zero for \z\ — > oo, V X x tends asymptotically (as — 1/|#|) towards the 
constant (positive) contribution of Eq.(H). This results in larger intersubband spacings 
using XX theory for the same total electron density. Note that the minimum of V xx (z,a) 
is not at the maximum of n(z, a) as in V£ D {z, a) and also note that V xx (z,a) is large 
where n (z, a) is negligible. Both features are consequences of the non-local dependence of 
V xx (z, a) on n (z, a) , displayed explicitly in Eq.(6). 

Significant qualitative differences between LDA and XX theory appear each time a new 
subband is occupied. For instance in Eq.(|^), at zero temperature, a finite term proportional 
to 5Ex/5k l p appears disco ntinuously when k l p — > + . Moreover, in Eqs.(|) and (^|), 
when /i — + el + + , a finite contribution (— \(fi(z, z')\ 2 /2nA) appears in the second term. 
Therefore, the inverse \a l { z i z> ) a ^ so changes discontinuously when k l p — > + . Although 
the discontinuities in SEx/dk}? and Xu 1 have opposite effects, they do not cancel each other 
and Xa 1 dominates [Q. Let us consider what happens when /z crosses a subband energy 
e a n but one neglects the self-consistent adjustments of the charge. In that case, because 
of the discontinuities introduced by 5Ex/Sk F y and in x" 1 , the XX potential must change 
discontinuously: 

V xx (z, a, fi - el + 0+) = V xx (z, a^^e^- 0+) + AV^iz, a). (7) 

Therefore, the wave functions, the total electron density, and the total energy cannot be 
continuous functions of /x and the discontinuity AV XX ' n signals a first-order transition of 
the 2DEG. An interesting question is what is the effect of AVx X (z, a) on the intersubband 
spacing E i. Case I: if one only considers the discontinuities in SEx/Sk}?, E 01 (fi — > el + 
+ ) < E i(fi — > el — + ) and \x lies above the first-excited subband and a finite amount of 
charge is abruptly transferred from the ground to the first-excited subband. Case II: if one 
adds the discontinuities in x" 1 , E i(fi — > e\ + + ) > E i(fi — >• e\ — + ) and \i remains below 
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the bottom of the first-excited subband and a self-consistent solution of the system is not 
possible (under our assumption of translational symmetry). In words, if the system attempts 
the occupation of the first-excited subband, this increases the intersubband spacing, which 
in turn empties the excited subband in the next iteration, which in turn produces a lower 
intersubband spacing, and so on. Finally, if E i(fi — > e\ + + ) = E 01 (fi — ■> e\ — + ) (Case 
III) a smooth occupation of the first-excited subband would be possible. Provided that self- 
consistency is frozen, the only possibility in the LDA is Case III, because the LDA exchange 
potential is a continuous function of the density. Self-consistent effects, in some conditions, 
can generate first-order transitions within LDA ||. But in XX theory, in contrast, a second 
order transition will be an accident and, in general, first-order transitions are expected when 
H crosses a subband energy at zero temperature. [|1^ 



In Fig.^ we plot the intersubband spacing E m as a function of the total 2D electron 
density N/A = (iVj + N±)/A with different methods. The dotted line corresponds to the 
usual LDA, the continuous line to XX theory at zero temperature, and the dashed line was 
obtained with an approximated [|TTJ finite temperature XX theory (T = 10 K in Eq.(|J)). 
The straight line is the evolution of fi-e%= (7ih 2 /m*)N/A in an hypothetical system with 
a single band occupied. All methods give similar qualitative results when E Q i is far from 
fi — Eq, though LDA calculations give smaller intersubband spacings than XX methods (as 
expected from Fig.||]). However, significant qualitative differences appear when E m « H~Eq. 
The LDA calculation gives a continuous curve but there is a discontinuity in the derivative 
of Eqi when the first-excited subband is occupied (which implies a second order transition). 
The zero temperature calculations with XX theory give a remarkable result: there is a 
window of densities where a self-consistent solution is not achieved because it corresponds 
to the Case II discussed above. Our interpretation of this window, where it is not possible to 
find a self-consistent solution under the assumption of translational invariance, is that this 
symmetry must be broken in the ground state. Thus in this density region a new broken 
symmetry phase might exist at low temperatures. 

Consideration of finite temperatures in Eq.(||) introduces smoothly the discontinuity in 
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Xa- 1 (z, z') when \fj, — e°\ « fc^T. This allows us to achieve a self-consistent solution for 
all densities (although the numerical convergence becomes unstable in some cases). As 
the density increases, the occupation of the first-excited subband cannot be avoided and 
5Ex/5k l p generates a sudden drop in Eq\. This new solution is also stable at lower den- 
sities. Accordingly, there is a range of densities where it is possible to find two solutions. 
A first-order transition occurs when the free energies of the two solutions cross; detailed 
temperature-dependent calculations will be reported elsewhere. From Fig|| we can estimate 
the drop on Eqi considering the energies in the density range where the two solutions exist: 
Ai?oi = 2.4 ± 0.4 meV, which is in very good agrement with the measurement of 3.5 meV 
reported by Goni et al. ||. 

In summary we have extended the KS-DFT for 2DEG beyond the state of the art method 
based on the LDA. Our theory allows to obtain an exact exchange potential without ap- 
proximations up to numerical precision, while correlations are considered in the LDA level. 
The theory satisfies known limits and theorems. The theory predicts first-order transitions 
every time a new subband is occupied and suggest that for some systems at zero tempera- 
ture the translational symmetry might be spontaneously broken. Finite temperature effects 
are included approximately. We calculated the self-consistent solutions for a realistic sys- 
tem and obtained phase transitions which are in good quantitative agreement with recent 
experiments. 
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FIGURES 

FIG. 1. (a) Schematic representation of an asymmetric doped quantum well nanostructure. (b) 
Self-consistent subband structure and self-consistent potential for the zero-bias situation. 

FIG. 2. Comparison of the exact exchange potential with the LDA potential when only one 
subband is occupied. The full upper line is the 3D density n(z) obtained with XX theory (corre- 
sponding to a 2D density of 0.68 x 10 12 /cm 2 ). XX and LDA exchange potentials are compared for 
this same density. Energies are measured from the chemical potential \i. 

FIG. 3. Intersubband spacing Eq\ as a function of the 2D density N/A. Full lines, zero tem- 
perature exact exchange; dashed lines, finite temperature exact exchange; dotted line, LDA. 
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